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Summary 


This  report  summarizes  the  results  of  the  research  and  development  effort  under  the  auspice 
of  AFRL  Award  |jFA8750-l  1-1-0040.  The  primary  objective  of  the  project  is  to  study  and 
implement  MIMO  architecture  in  airborne  communications.  The  primary  challenges  include 
1)  the  lack  of  scatterers  in  airborne  communications  which  often  renders  the  channel  matrices 
to  be  rank  dehcient;  2)  fast  channel  variation  due  to  high  mobility  of  the  communication 
platforms. 

Toward  overcoming  these  challenges  and  in  close  collaboration  with  AFRL  researchers, 
we  have 

•  explored  the  feasibility  of  channel  tracking  in  a  dynamic  channel  environment; 

•  developed  a  GUI  system  that  allows  one  to  visualize  the  effective  capacity  under  various 
communication  scenarios; 

•  developed  a  D-BLAST  architecture  that  employs  orthogonal  spreading  to  counter  the 
effect  of  channel  rank  dehciency; 

•  implemented,  using  python  and  with  four  USRP  N210  [1]  nodes,  a  2  x  2  MIMO  system 
with  over  the  air  transmission. 

The  following  refereed  articles  have  been  published  under  this  effort: 

•  K.  Borle,  B.  Chen,  and  M.  J.  Cans,  Channel  tracking  for  D-BLAST  for  airborne 
platforms,  Proc.  Asilomar  Conference  on  Signals,  Systems,  and  Computers,  Monterey, 
CA,  Nov.  2011 

•  M.  J.  Cans,  K.  Borle,  B.  Chen,  T.  Freeland,  D.  McCarthy,  R.  Nelson,  D.  Overrocker 
and  P.  Oleski,  Enhancing  connectivity  of  unmanned  vehicles  through  MIMO  communi¬ 
cations,  Proc.  of  IEEE  Vehicular  Technology  Conference,  Las  Vegas,  NV,  Sept.  2013. 

In  addition,  various  software  packages  have  been  delivered  to  the  AFRL  collaborators, 
including 

•  Channel  tracker  simulator,  as  described  in  Chapter  3. 
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Channel  matrix  simulator.  This  is  to  take  the  Newport  measurement  data  and  con¬ 
struct  a  new  channel  matrix  according  to  the  relative  position  of  the  transmitter  and 
receiver  through  extrapolation. 

LDPC  encoder  and  decoder  source  code. 

C-|--f  code  for  the  MIMO  D-BLAST  transmitter  and  receiver  blocks. 
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Introduction 


A  point  to  point  airborne  MIMO  communication  link  can  be  defined  as  two  terminals  com¬ 
municating  over  a  wireless  channel  such  that  either  one  or  both  the  terminals  are  airborne 
and  both  employ  multiple  antennas.  On  the  other  hand,  the  conventionally  studied  MIMO 
wireless  communication  is  based  on  the  premise  that  the  communicating  multi-antenna  ter¬ 
minals  are  located  in  a  dense  urban  environment.  Two  important  factors  distinguish  such  an 
airborne  MIMO  channel  from  a  conventional  MIMO  channel:  Firstly,  the  channel  conditions 
change  very  rapidly  in  an  airborne  link  due  to  fast  relative  motion  between  the  airborne  ter¬ 
minals.  Secondly,  the  channel  offers  a  significantly  less  scattering  environment  as  compared 
to  a  dense  urban  environment.  As  such,  airborne  MIMO  communication  poses  a  somewhat 
different  set  of  problems  as  compared  to  the  conventional  MIMO  communication  problems. 

This  report  describes  our  efforts  in  addressing  the  above  mentioned  problems  of  airborne 
MIMO  communication.  To  be  precise  we  address  the  problems  of  channel  tracking  in  dy¬ 
namic  channel  conditions  and  communication  over  rank  deficient  channels.  In  addressing 
these  problems,  we  have  encountered  D-BLAST  transceiver  architecture  as  the  central  idea, 
upon  which  our  proposed  solutions  are  based.  Hence,  to  facilitate  experiments,  we  also  de¬ 
velop  a  GNU  Radio/USRP  based  D-BLAST  communication  system  employing  2  transmit 
and  2  receive  antennas.  Furthermore  we  have  developed  a  C-I--I-  library  that  provides  D- 
BLAST  encoding  and  decoding  functionality.  The  C-I--I-  library  is  not  coupled  to  the  GNU 
Radio  framework  and  hence  provides  more  flexibility  in  its  usage. 

In  the  following  we  provide  a  brief  description  of  the  contributions  followed  by  the  orga¬ 
nization  of  the  report. 

2.1  Channel  Tracking 

The  channel  in  airborne  communication  changes  continuously  at  a  rapid  rate.  To  facilitate 
coherent  symbol  detection  in  such  a  rapidly  changing  channel,  the  receiver  needs  to  update 
the  channel  state  information  at  a  much  more  frequent  rate.  The  traditional  use  of 
embedding  pilot  symbols  with  payload  data  may  not  be  easy  to  justify  if  channel 
variation  is  fast  and/or  data  rate  requirement  is  high.  We  consider  the  implementation 
of  D-BLAST  for  airborne  platforms  and  develop  channel  tracking  scheme  to  eliminate  or 
reduce  the  use  of  pilot  symbols.  In  a  normal  operation  mode,  channel  update  is 
achieved  dynamically  as 
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each  layer  of  the  D-BLAST,  encoded  using  an  LDPC  code,  is  decoded.  To  ensure  that 
the  transceiver  can  detect  outage  due  to  the  loss  of  channel  state,  an  adaptive  algorithm  is 
devised  utilizing  the  extremal  property  of  terminating  likelihood  ratio  of  an  LDPC  decoder. 


2.2  Variable  Rate  MIMO 

In  this  part  of  the  report  we  develop  a  variable  rate  MIMO  scheme,  based  on  D-BLAST 
transceiver  architecture,  to  overcome  MIMO  channels  that  are  rank  dehcient.  Real  channel 
measurements  are  conducted  and  the  analysis  supports  the  use  of  MIMO  communications  for 
such  applications  due  to  its  potential  throughput  advantage.  Unique  challenges  and  the  ways 
to  address  them  are  described  in  detail.  In  particular,  the  lack  of  scattering  and  the  blockage 
of  line  of  sight  may  lead  to  rank  dehcient  channel  matrices,  which  are  exacerbated  due  to 
the  absence  of  channel  state  information  at  the  transmitter.  A  variable  rate  MIMO  scheme 
is  thus  proposed  to  overcome  these  challenges  in  order  to  realize  the  promising  throughput 
gain  afforded  by  MIMO  communication. 


2.3  GNU  Radio  MIMO  Implementation  Using  USRP 
N210 

An  indispensable  part  of  any  scientihc  inquiry  is  its  experimental  component.  The  afore¬ 
mentioned  Variable  Rate  MIMO  scheme  shows  promising  results  for  the  rank  dehcient  chan¬ 
nels  obtained  through  real  measurement  data.  The  results,  as  encouraging  as  they  appear, 
nonetheless  must  be  validated  through  experimental  means.  Therefore,  to  conduct  real-time 
wireless  experiments  we  need  a  platform  that  facilitates  efficient  and  fast  prototyping  of  a 
testbed.  Software  dehned  radio,  which  provides  us  with  these  desirable  traits,  appears  as  a 
very  attractive  option  for  developing  the  necessary  components  required  for  the  experimen¬ 
tation. 

The  GNU  Radio  software  development  toolkit  combined  with  USRPs  provide  an  effi¬ 
cient  and  rapidly  employable  software  dehned  radio  system.  Using  the  GNU  Radio/USRP 
platform,  we  have  developed  a  2  x  2  D-BLAST  transceiver.  We  have  also  created  a  G-I--I- 
library,  developed  separately  from  the  GNU  Radio  based  testing  framework  but  able  to  be 
used  within  the  framework. 


2.4  Organization 

The  rest  of  the  report  is  organized  as  follows. 

•  Ghapter  3  develops  and  analyzes  the  channel  tracking  scheme  for  a  D-BLAST  based 
airborne  communication  system. 

•  Ghapter  4  develops  a  variable  rate  MIMO  scheme,  based  on  D-BLAST  transceiver 
architecture,  to  overcome  MIMO  channels  that  are  rank  dehcient. 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 

4 


Chapter  5  provides  a  brief  overview  of  the  C++  based  software  library  which  is  devel¬ 
oped  to  facilitate  experimentation  of  MIMO  communication. 

Chapter  6  concludes  the  report. 
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Channel  Tracking 


3.1  Introduction 

For  a  communication  system  in  which  either  one  or  both  the  stations  are  in  motion,  the 
channel  is  continuously  changing.  At  the  receiver  end,  for  coherent  detection  of  symbols, 
this  translates  to  constant  updating  of  channel  state  information.  The  conventional  method 
is  to  multiplex  pilots  and  data  symbols.  For  the  case  where  there  is  rapid  variation  in 
channel,  this  pilot  overhead  may  not  be  affordable.  In  order  to  tackle  this  scenario  we  try 
to  investigate  the  case  in  which  we  use  previously  detected  symbols  to  estimate  channel  for 
next  block. 

Since,  the  channel  varies  temporally,  simply  using  the  previous  block  of  decoded  symbols 
may  not  be  optimal.  We  try  to  solve  this  problem  by  hrst  hnding  the  optimal  length  of 
symbols  needed  to  minimize  the  mean  square  error  of  estimated  channel  for  a  SISO  channel 
with  rayleigh  fading.  Then  we  generalize  this  to  MIMO  case. 

Once  we  know  how  to  track  the  channel  efficiently,  the  question  arises,  how  to  detect 
outage  at  the  receiver.  We  narrow  down  our  search  for  an  outage  detection  scheme  to  systems 
which  use  LDPC  codes.  We  hnd  through  simulations  that  the  average  mean  magnitude  of 
the  log  liklihood  ratios  remain  small  when  the  codeword  is  docoded  incorrectly,  while  in  the 
otherwise  case,  the  mean  magnitude  has  a  high  value. 

Section  3.2.1  gives  the  analysis  for  optimal  length  to  minimize  the  mean  squared  error  in 
a  Rayleigh  fading  channel.  We  then  give  a  brief  description  of  the  D-BLAST  architecture  in 
Section  3.2.2.  Section  3.2.3  describes  the  channel  tracking  algorithm  for  D-BLAST,  followed 
by  simulation  results.  Section  3.4  briefly  describes  the  channel  tracking  simulator,  which  we 
have  developed. 

3.2  Methods,  Assumptions  and  Procedures 

3.2.1  Channel  Estimation  -  Optimal  Length 

Single  Input  Single  Output  Case 

Consider  a  SISO  channel,  yi  =  hiXi  +  Zi,  where  Xi  is  the  transmitted  symbol,  Zi  is  the  additive 
noise  and  hi  is  the  fading  coefficient  in  the  time  slot  [2].  If  the  transmitter  transmits  A^-|- 1 
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symbols  in  iV  +  1  consecutive  time  slots, 
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Write  the  diagonal  matrix  of  fading  coefficients  as  h^-k  =  +  A^,  for  /c  =  1,2, . . .  iV 

[3],  eq.  (1)  becomes 
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where  I^r  is  an  iV  x  iV  identity  matrix  and  x  =  [x^  Multiplying  both  sides  by 

we  get, 

N  N 

X^y  =  /lArX^X  +  ^  ^kXkxl  +  ^  xlzk.  (3) 

fc=l  fc=0 

For  simplicity,  let  XkX^.  =  1.  The  linear  least  squares  estimate  is  given  by. 
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The  error  in  estimating  Hn  thus  arises  from  two  factors:  One  due  to  temporal  channel 
variation  and  the  other  due  to  additive  noise  component.  Our  goal  is  to  hud  the  optimum 
training  length  that  provides  the  best  trade-off  between  these  two  error  types.  This  is 
equivalent  to  minimizing  the  average  of  |/iAr  —  /lArp.  To  achieve  this  we  hrst  try  to  hnd  a 
tractable  form  of  the  MSE, 


/lATp 


=  S  (Jin  —  (J^N  —  h]\f 
1 


{N  +  iy 


N  \  /  ^  ^ 
k=l  /  \k=l  / 


+ 


N  +  V 


(5) 

(6) 


where  is  the  noise  variance.  The  expectation  of  the  cross-product  terms  is  0  because 
we  assume  that  the  noise  is  zero  mean  and  is  independent  of  channel  coefficients.  To  evaluate 
^l/iAT  —  /iTvP,  we  hrst  evaluate  the  expectation  in  the  hrst  term  in  (6).  We  use  the  approach 
similar  to  [4] .  We  assume  the  channel  as  frequency  hat  fading  with  a  dihuse  scattering  channel 
model,  which  lends  itself  mathematical  tractability.  The  following  theorem  summarizes  our 
result. 
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Theorem  1.  Define  =  h{{N  —  k)Ts)  —  h{NTs).  h(t)  is  the  channel  modeled  as  follows, 


h(t)  =  f;  (7) 

k=l 

where  S  is  the  number  of  scatterers,  fik  and  fk  are  the  complex  amplitude  and  Doppler 
frequency  of  the  multipath  respectively  and  Tg  is  the  symbol  period,  with  the  following 
assumptions. 

1.  13 k  are  independent  zero  mean  random  variables,  and  normalized  such  that  I]f=i  l/SfcP  = 

1  . 

2.  fk  are  i.i.d.  distributed  uniformly  on  the  interval  [— /d,/d]- 

3.  fik  and  fk  are  independent  of  each  other. 


Then,  (See  Appendix  A) 


=  N{N  +  1) 


smijifoNTs)  sin(7r/D(7V  +  1)T^) 

nfoTg  sm{7r fnTs) 


Therefore,  putting  (8)  in  (6),  we  get 


(8) 


<o\hj\f  —  /ijv| 


N 

iV  +  1  ^  iV+  1 

sm^nfpNTg)  sin(7r/g(7V  +  1)TJ 
(A^  +  1)%/dT^  sm{7r  fpTg) 


(9) 


The  above  expression  gives  the  MSE  as  a  function  of  maximum  Doppler  spread,  symbol 
period  and  number  of  training  symbols.  For  a  given  Doppler  spread  and  symbol  time,  we  can 
find  the  optimal  training  length  numerically.  We  consider  the  MIMO  case  in  the  following 
section. 

A  plot  of  MSE  against  the  fraction  of  code  block  size  used  for  estimation  is  given  in 
Fig.  1  for  fpTg  =  1.25  x  10“"^,  where  fp  is  the  maximum  Doppler  spread  and  Tg  is  the 
symbol  period. 


Multiple  Input  Multiple  Output  Case 

Consider  a  MIMO  system  with  t  transmit  antennas  and  r  receive  antennas.  We  assume  the 
following  channel  model,  y*  =  HjXj  +  Zj,  where  x*  is  the  f  x  1  transmitted  vector,  is  the 
r  X  1  received  vector,  Hj  is  the  r  x  t  channel  matrix  and  Zj  is  the  r  x  1  noise  vector  at  time 
index  i.  We  assume  the  noise  to  be  complex  Gaussian  with  zero  mean  and  covariance  matrix 
equal  to  the  identity  matrix. 

Let  the  transmitter  send  iV  + 1  symbol  vectors  in  iV  + 1  time  slots.  The  iV  + 1  transmitted 
symbols  from  one  antenna  are  orthogonal  to  that  of  other  transmit  antennas. 
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- Closed  Form 

- Monte  Carlo 


0  0.2  0.4  0.6  0.8  1 

Fraction  of  Code  Block  (Block  Size  =  200  Symbols) 

Figure  1:  Behavior  of  MSE  of  SISO  channel  for  foTg  =  1.25  x  10“^  with  block  length  of  200 
symbols 


Yn  =  HnXn  +  zn, 

YN-I  =  Hn_iXn_i  +  Zn_i 
=  HnXn_i  +  AiXn_i 


ZN-1, 


(10) 


(11) 


Yo  =  HqXo  +  zo 

=  HnXo  +  AnXo  +  Zq. 

The  MSE  in  estimating  Hat  can  be  given  as  follows, 


(12) 


^l|H 


N 


Hiv|||  =  ^E  V 


(ti) 


h 


(iJ) 

N 


=  E^ 


ihj)  _ 

Hn  iI'N 


=  rt<g 


h 


(1,1) 

N 


h 


(1,1) 

N 


is  the  channel  coefficient  between  the  receiver  antenna  and  the  transmitter 
antenna  of  Hn.  The  third  equality  is  because  all  the  elements  are  assumed  i.i.d.  The 
expectation  in  (15)  can  be  reduced  to  that  in  expression  (9). 

Therefore,  we  can  write  the  MSE  as 

(  (^1  N 

=  rt  - 5 - ^ - 

I  iV  +  1  iV  +  1 


(13) 

(14) 

(15) 


sm{7r fpNTs)  sm{7rf d{N  +  1)TJ ' 
(TV  +  1^71  fpTs  sm{7r  fpTs) 


(16) 
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Fraction  of  Block  (Block  Size  =  324  Symbols) 

Figure  2:  Behavior  of  MSE  of  a  4  x  4  channel  for  /dT*  =  7.72  x  10“^  with  block  length  of 
324  symbols. 

From  Fig.  2  we  see  that  for  about  30  percent  of  a  code  block  of  length  324  symbols,  lowest 
MSE  is  achieved  for  fnTg  =  7.72  x  10“^.  We  can  use  this  as  a  guideline  to  set  the  length  of 
symbols  in  our  tracking  algorithm.  The  tracking  algorithm  is  explained  in  Section  3.2.3. 

Fig.  3  shows  the  MSE  for  our  channel  tracking  algorithm.  For  a  setup  similar  to  above, 
the  lowest  MSE  for  channel  estimation  is  achieved  when  only  about  30  percent  of  code  block 
is  used.  The  next  section  gives  a  brief  explanation  of  the  D-BLAST  architecture  and  then 
the  channel  tracking  algorithm. 


Number  of  Code  Blocks 


Figure  3:  4  x  4  D-BLAST  tracking  MSE  w.r.t  sub-blocks  for  foTs  =  7.72  x  10 

3.2.2  D-BLAST  Architecture 

D-BLAST  stands  for  Diagonal  Bell  Labs  Space  Time  Architecture.  It  is  a  transceiver  archi¬ 
tecture  for  MIMO  systems  [5].  For  an  r  x  t  (r  receive  antennas  and  t  transmit  antennas) 
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MIMO  system,  the  input  data  stream  is  demultiplexed  into  t  streams.  We  call  each  de¬ 
multiplexed  stream  as  a  sub  stream.  The  t  sub  streams  are  then  independently  coded  and 
modulated.  These  t  sub  streams  are  then  transmitted  over  the  t  transmit  antennas.  However, 
the  association  between  the  sub  streams  and  the  antennas  is  changed  in  a  periodic  manner. 
Each  sub  stream  is  rotated  through  the  t  transmit  antennas  during  one  code  block  duration. 
This  results  in  diagonal  layering  of  the  sub  streams  across  the  antennas  and  time.  At  the 
receiver  end,  decoding  happens  one  layer  at  a  time.  Because  of  the  diagonal  structure,  each 
sub  stream  sees  channels  from  all  the  transmit  antennas  and  hence  has  improved  diversity 
compared  with  V-BLAST  [6].  This  results  in  spreading  out  the  errors,  in  the  event  that 
one  of  the  links  has  a  bad  channel  condition,  thereby  minimizing  the  probability  of  out¬ 
age.  The  diagonal  structure  renders  D-BLAST  as  an  outage  optimal  scheme  in  slow  fading 
environments  [2]. 

3.2.3  Channel  Tracking 

For  the  sake  of  illustration  we  describe  our  scheme  for  a  2  x  2  MIMO  system.  It  can  be 
extended  to  any  q  x  q  MIMO  system  in  a  straightforward  manner.  The  symbol  space-time 
diagram  is  shown  below.  The  stream  of  symbols  Xjj  denote  that  they  belong  to  stream 
i  and  are  transmitted  during  the  time  block,  tk  and  ak  denote  the  time  block  and 
antenna  number  respectively,  over  which  the  indexed  symbols  are  transmitted.  We  use  the 
jth  block  abstraction  to  indicate  the  time  duration  during  which  the  transmitter  keeps 
the  association  between  the  streams  and  the  antennas  constant,  p  indicates  pilot  symbols. 
li  denotes  layer.  For  example:  xiiXi2  is  layer,  X21X22  is  2”'^  layer  and  so  on. 


to 

tl 

t2 

to 

ai 

Pai 

Xll 

X2I 

Xl3 

0,2 

Paa 

Pa2 

X12 

X22 

Table  1:  Symbol  Space-Time  Diagram 

Let  Yj  and  X*  be  the  received  and  transmitted  set  of  vector  in  time  block  i.  We  assume 
that  we  already  have  the  initial  channel  estimate  Hq,  found  during  time  block  tossing  the 
pilots  Pai  and  p^j  as  shown  in  previous  section.  We  also  assume  that  the  pilot  p^j  has  been 
subtracted  out  from  the  received  vector  during  tj. 

In  this  scheme,  we  use  the  initial  estimate  Hq  to  detect  xn  and  X12.  The  receiver  then 
decodes,  re-encodes  and  performs  baseband  modulation  to  reconstruct  xn  and  X12.  xnalong 
with  pa2  are  used  to  hnd  a  new  estimate  Hi  of  the  channel.  The  tensor  product  of  X12  and  its 
corresponding  channel  vector  from  Hi  is  then  subtracted  from  Y2,  where  Y2  is  the  channel 
output  during  the  time  block  t2.  This  new  channel  estimate  Hi  is  then  used  to  detect  the 
next  layer  i.e.,  X21  and  X22,  and  the  process  is  repeated. 
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3.3  Results  and  Discussion 


3.3.1  Channel  Tracking 

For  simulation  purposes  we  consider  a  4  x  4  MIMO  system  with  a  Rayleigh  fading  channel. 
For  error  control  coding  we  use  an  LDPC  (Low  Density  Parity  Check)  code  of  block  length 
1944.  The  bit  streams  are  then  modulated  by  4-QAM  modulators.  We  show  simulation 
results  for  2  cases  of  channel  variations  and  2  cases  of  code  rates. 


- f„T  =15.43..  10'* 

D  s 

- UT  =7.72:.:10'' 

D  s 

J 

/' 

} 

1 
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I 

II 
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1 
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■’I _ 1 _ I _ I 
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Number  of  Code  Blocks 


Figure  4:  Behavior  of  channel  estimation  mean  square  error  with  increasing  number  of  code 
blocks 
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Figure  5:  Behavior  of  the  mean  magnitude  of  LLR  of  decoded  LDPC  code  with  increasing 
number  of  code  blocks 


We  can  see  that  for  a  faster  varying  channel,  the  algorithm  loses  track  very  early.  The 
sudden  jumps  in  MSE  in  Figure  4  indicate  loss  of  channel  tracking.  They  correspond  to  a 
sudden  change  in  the  mean  magnitude  of  LLRs  of  decoded  LDPC  codes  in  Figure  5.  We  can 
use  these  sudden  changes  to  detect  outage  at  the  receiver. 
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Figure  6:  Behavior  of  Average  Mean  Magnitude  of  Log  Likelihood  Ratios  for  an  LDPC  code 
with  N  =  648  and  R  =  ^ 

3.3.2  Outage  Detection 

During  the  tracking  process,  decoding  errors  may  occur.  The  incorrectly  decoded  codeword 
will  lead  to  an  inaccurate  channel  estimate.  In  order  for  the  receiver  to  detect  this,  and 
thereby  tell  the  transmitter  to  reinitialize  the  transmission,  we  use  the  mean  magnitude 
of  LLRs  of  LDPC  codes  [7].  We  found  throngh  simnlations  (see  Fig.  6)  that  the  mean 
magnitnde  of  LLRs  of  correctly  decoded  LDPC  codes  increases  with  the  nnmber  of  message 
passing  iterations.  In  comparison,  for  incorrectly  decoded  codewords,  the  mean  magnitnde 
remains  constant  at  some  small  valne.  We  can  use  this  property  to  detect  the  loss  in  tracking. 
Loss  in  tracking  resnlts  in  incorrectly  decoded  codewords  or  vice  versa,  thereby  resnlting  in 
smaller  mean  magnitnde  of  the  LLRs.  Using  a  threshold  detector  we  can  detect  this  loss  in 
tracking  and  thus  detect  outage. 

3.4  Channel  Tracking  Simulator 

In  the  previous  sections  we  stndied  channel  tracking  for  a  D-BLAST  commnnication  system 
nsing  payload  data.  In  this  section  we  provide  an  overview  of  the  tool  we  developed  to 
simulate  such  channel  tracking.  Using  this  tool  we  can  visnalize  the  performance  of  the 
channel  tracking  algorithm  for  different  set  of  parameters  with  ease.  The  tool  was  developed 
in  MATLAB  [8]. 

3.4.1  Component  Description 

There  is  only  one  major  component  to  the  GUI  -  the  main  window  as  shown  in  7.  R  is  used 
to  input/select  different  parameters,  which  are  then  passed  on  the  appropriate  fnnctions  for 
simnlation.  We  give  a  brief  description  of  the  components  of  the  main  window  below. 

•  ^  of  Transmit  Antennas  is  use  to  set  the  number  of  transmit  as  well  as  receive 
antennas  we  assnme  in  ont  calculations  that  they  are  equal  in  nnmber.  Inputting  a 
non-integer  valne  will  cause  an  error. 
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Figure  7:  Channel  Tracking  Simulator 


•  7^  of  Receive  Antennas  is  not  editable.  Its  value  is  set  equal  to  that  of  Transmit 
Antennas. 

•  Carrier  Frequency  is  used  to  set  the  carrier  frequency  in  Hertz. 

•  Bandwidth  is  used  to  set  the  communication  bandwidth  in  Hertz. 

•  Modulation  Order  is  used  to  select  either  BPSK  (2)  or  QPSK  (4)  mode  of  modula¬ 
tion. 

•  No  Channel  Coding  sets  the  simulation  to  not  to  use  channel  coding.  If  this  item 
is  selected  the  following  item  LDPC  cannot  be  selected. 

•  LDPC  sets  the  simulation  to  use  LDPC  codes  during  simulation.  There  are  in  all  6 
combinations  while  using  LDPC  codes.  Code  blocks  of  lengths  648,  1296  and  1944. 
While  code  rates  options  are  3/4  and  1/2.  The  parity  check  matrices  for  these  codes 
are  obtained  from  802.1  In  standard  [9]. 

•  Don’t  Track,  if  selected  doesn’t  perform  any  tracking.  It  is  suitable  only  if  we  keep 
the  channel  constant. 

•  Track  enables  the  selected  scheme  for  channel  estimation  or  tracking,  depending  on 
the  scheme.  The  scheme  descriptions  are  given  below. 

—  Scheme  1  -  the  tracking  scheme  given  in  section  3.2.3. 

—  Scheme  2  -  pilots  multiplexed  with  data  in  the  initial  t  —  1  sub  blocks. 
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—  Scheme  3  -  same  as  scheme  1,  except  that  in  this  case  we  can  set  how  much 
fraction  of  code-block  to  be  used  for  channel  tracking,  whereas  in  scheme  it  was 
hxed  to  one  sub  block. 

—  Scheme  4  -  a  scheme,  wherein  we  distribute  the  pilots  across  the  payload. 

—  Fraction  of  Code  Block  -  this  value  has  different  interpretations.  For  Scheme 
1  &  2  it  is  of  no  use.  For  Scheme  3,  it  means  the  fraction  of  code  block  to  be  used 
for  channel  estimation.  In  this  case  we  can  set  it  as  an  array  of  different  values. 
For  Scheme  4,  it  indicates  the  fraction  of  code  block  between  two  pilot  symbol 
blocks.  The  value  should  not  be  an  array  in  this  case. 

•  Maximum  Doppler  Shift.  This  edit  box  is  not  functional. 

•  Gamma  is  used  to  set  the  maximum  Doppler  shift.  Gamma  is  dehne  as,  7  = 

time  to  send  one  code  block _  gamma  values  to  Calculate  maximum  Doppler  shift, 

which  are  then  used  to  generate  the  channel  coefficients. 

•  Transmit  Block  Length  in  ^  of  Codeword  Lengths  is  self-explanatory. 

•  Pilot  Length  is  used  to  set  the  length  of  pilot  block  in  number  of  symbol  periods. 

•  SNR  is  used  to  set  the  average  signal  to  noise  ratio  at  each  receiver  antenna  in  decibels. 

•  ^  of  samples  to  be  averaged  are  the  number  of  Monte  Carlo  runs  that  we  average 
to  get  the  sample  statistics. 


3.5  Conclusion 

In  this  chapter  we  studied  channel  tracking  for  a  D-BLAST  communication  system  using 
payload  data.  Since  the  tracking  involves  using  previously  detected  data  symbols,  a  closed- 
form  expression  for  channel  MSE  as  a  function  of  the  number  of  training  symbols,  maximum 
Doppler  spread  and  symbol  duration  was  derived  for  a  SISO  link  in  a  diffuse  scattering 
environment.  This  expression  was  then  extended  to  MIMO  links.  We  use  this  analysis  to 
hnd  the  optimum  length  of  previously  detected  symbols  for  channel  tracking.  In  the  event 
the  receiver  loses  track  of  the  channel,  we  show,  using  simulations,  that  using  the  mean 
magnitude  of  log  likelihood  ratios  of  LDPC  codes  can  detect  the  event  of  a  decoding  error, 
which  in  turn  can  be  used  to  detect  outage. 

We  have  also  provided  a  MATLAB  based  channel  tracking  simulator.  It  can  be  used  to 
easily  simulate  the  developed  channel  tracking  and  outage  detection  schemes. 
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Figure  8:  Scheme  1  Flowchart.  U  =  [si^jS2,i+i . . .  Hj 
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Variable  Rate  MIMO 


4.1  Introduction 

The  use  of  autonomous/unmanned  vehicles  for  various  civilian  and  military  applications 
has  become  increasingly  prevalent.  The  absence  of  human  pilots  on  board  those  unmanned 
systems,  however,  make  it  a  challenging  task  to  accomplish  various  intended  missions.  For 
example,  these  unmanned  systems  are  often  remotely  piloted  and  thus  having  a  reliable  and 
robust  communication  link  between  the  aerial  systems  and  the  ground  control  unit  is  impera¬ 
tive.  Even  for  systems  that  are  semi-autonomous,  having  a  high  throughput  communication 
link  is  often  essential  to  accomplishing  any  intended  mission.  Many  remotely  piloted  air¬ 
crafts  (RPAs)  are  used  for  surveillance  applications  and  the  need  to  stream  surveillance  data, 
including  real  time  video  data  requires  a  highly  reliable  and  high-throughput  communication 
link  from  the  RPAs  to  ground  units. 

Many  legacy  communication  links  (e.g..  Link  16  with  a  data  rate  not  more  than  16kb/s) 
operate  at  a  data  rate  that  becomes  highly  inadequate  for  applications  where  video  streaming 
from  aerial  systems  to  ground  is  needed.  Merely  scaling  up  the  power/bandwidth  is  both 
limited  by  resource  and  policy  constraints  as  well  as  the  fundamental  theoretical  limits 
dictated  by  the  Shannon  theory.  A  promising  technology  is  the  use  of  multiple  antenna 
communication  systems  [5,10,11].  The  so-called  multiple-input  multiple-output  (MIMO) 
communication  scales  up  data  rate  linearly  as  the  number  of  antennas  increase  and  thus 
provides  great  potential  for  improving  the  throughput  of  air  to  ground  communication.  This 
helps  enable  many  envisioned  applications  that  may  otherwise  be  infeasible. 

MIMO  communications,  however,  are  traditionally  designed  for  the  so-called  scattering 
environment  [2]  where  independent  channel  variations  between  different  transmit /receive 
antenna  pairs  are  exploited.  For  airborne  platforms,  however,  there  has  been  a  debate  about 
the  feasibility  of  MIMO  communications  because  of  the  lack  of  scattering.  However,  for 
certain  communication  ranges,  the  large  aperture  that  an  aircraft  affords  makes  MIMO  an 
appealing  choice  of  communication  that  can  attain  signihcantly  higher  throughput  given  a 
hxed  power/bandwidth  budget  compared  with  single  antenna  systems  even  in  the  absence 
of  any  scatterers  [12]. 

This  chapter  describes  an  ongoing  research  and  development  effort  that  uses  MIMO 
communications  to  enable  robust  and  high  capacity  connectivity  between  RPAs  and  ground 
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terminals.  While  the  large  aperture  may  compensate  for  the  lack  of  the  scattering,  two 
unique  challenges  still  exist  for  airborne  MIMO  communications.  The  large  aperture  is  only 
attained  when  antennas  are  placed  strategically  apart  on  a  RPA.  In  the  absence  of  scattering, 
i.e.,  when  communications  are  limited  by  line-of-sight  channels,  the  fact  that  some  antenna 
elements  on  a  RPA  may  be  completely  out  of  sight  from  its  communicating  party  may 
render  the  channel  matrix  ill-conditioned.  This  is  further  complicated  by  the  high  mobility 
and  maneuverability  of  the  RPA  which  make  it  infeasible  to  have  complete  channel  state 
information  (CSI)  at  the  transmitter.  To  address  these  challenges,  we  propose  a  variable  rate 
MIMO  communication  scheme  that  combines  the  D-BLAST  architecture  with  per  antenna 
spreading  to  harvest  the  maximum  possible  throughput  gain  allowed  by  the  channel. 

The  chapter  is  organized  as  follows.  Section  4.2.1  describes  the  channel  measurement 
apparatus.  The  measurement  data  provides  guidance  on  the  potential  throughput  gain  for 
the  particular  application  of  interest  as  well  as  challenges  in  realizing  the  throughput  gain 
to  address  these  challenges.  Section  4.2.2  introduces  a  modihed  D-BLAST  architecture  to 
address  these  challenges.  Section  4.3  describes  simulation  results  using  the  real  measurement 
channels  to  show  the  potential  improvement  using  the  proposed  variable  MIMO  scheme. 

We  use  the  following  notations  throughout  the  chapter:  *  denotes  complex  conjugate, 
CN  stands  for  complex  Gaussian,  ^  is  transpose,  is  conjugate  transpose,  0  is  a  zero  vector 
and  \k  is  &  k  X  k  identity  matrix.  Small  letters  denote  scalars,  bold  small  letters  denote 
column  vectors  and  bold  capital  letters  denote  matrices. 

4.2  Methods,  Assumptions  and  Procedures 

4.2.1  MIMO  Channel  State  Matrix  Measurements 

Measurement  Apparatus 

The  process  of  measuring  the  channel  state  matrix  at  the  Newport  NY  radio  range  required 
measuring  the  complex  transmission  coefficient  from  each  antenna  element  on  the  ground 
array  to  each  element  on  the  unmanned  combat  air  vehicle  (UCAV)  model  of  a  RPA.  The 
UCAV  is  shown  in  Fig.  9. 


Figure  9:  UCAV  with  square  patch  antennas. 


The  UCAV  model  was  coated  with  Electrodag  to  shield  the  cable  networks  placed  inside 
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the  model  that  connect  the  patch  antennas  to  the  switch  which  feeds  the  RF  (Radio  Fre¬ 
quency)  receiver.  Although  this  shielding  is  not  necessary  in  an  operational  RPA,  it  helps,  in 
analyzing  the  measurements,  to  restrict  the  received  signals  to  the  patch  antennas  only.  Four 
patches  were  installed  for  the  uplink  band  (at  carrier  frequency  5.1  GHz  with  a  bandwidth 
of  100  MHz)  and  four  patches  for  the  downlink  band  (at  carrier  frequency  of  5.8GHz  with 
a  bandwidth  of  100  MHz).  Each  patch  has  two  ports  which  provide  perpendicular  linear 
polarization  with  about  30  dB  cross  polarization  discrimination  (isolation  between  ports). 
Thus  8  antenna  ports  are  provided  for  the  uplink  and  for  the  downlink,  respectively.  For  easy 
reference  these  ports  are  uniquely  labeled.  For  example,  RU12  represents  the  second  port  on 
the  hrst  uplink  patch  on  the  upper  side  of  the  UGAV.  The  antenna  patterns  were  measured 
on  the  Newport  Range.  Two  such  power  patterns  are  shown  in  Fig.  10(a)  and  10(b). 


((a))  Antenna  RUll  ((b))  Antenna  RU31 


Figure  10:  Antenna  Pattern 

The  power  pattern  for  antenna  RUll  gives  pretty  good  coverage  above  the  UGAV  because 
it  is  on  the  upper  surface  of  the  UGAV.  Antenna  RL31  is  on  the  bottom  surface  of  the  UGAV 
wing  and  has  good  coverage  in  the  hemisphere  below  the  UGAV  but  is  weak  above  the  UGAV. 
These  patterns  emphasize  the  need  for  the  D-BLAST  (Diagonal  Bell  Labs  Layered  Space 
Time)  [5]  form  of  MIMO  in  order  to  provide  full  coverage  for  all  data  streams. 

The  ground  array  utilizes  coax-to- waveguide  adapters  as  antenna  elements.  They  provide 
wide  angular  coverage  (110"  E-Plane  and  80°  H-Plane)  and  wide  bandwidth  4.9  GHz  to 
7.05GHz.  The  array  is  mounted  with  square  plates  as  shown  in  Fig.  11.  The  square  plates 
allow  a  choice  of  vertical  or  horizontal  polarizations.  Typically  measurements  are  made  with 
all  vertically  polarized  antenna  elements  or  with  polarizations  alternated  between  vertical 
and  horizontal  polarization  to  see  the  effect  of  polarization  mixing  on  MIMO  capacity.  The 
elements  can  be  at  various  spacing  by  bolting  the  plates  to  different  positions  along  the 
mounting  rails.  The  spacing  shown  in  Fig.  11  provide  an  80  inch  wide  array.  Another  test 
used  a  30  foot  wide  array.  The  array  is  mounted  on  an  expanding  tower  which  allows  array 
heights  of  2  feet  to  40  feet  above  ground  level. 

The  UGAV  is  mounted  on  a  positioning  tower  at  a  height  of  70  feet.  The  position  of  the 
UGAV  can  be  rotated  around  a  vertical  axis  (relative  to  ground)  for  360  degrees  of  azimuth. 
It  can  also  be  rotated  in  pitch  and  roll  over  90.  The  UGAV  support  is  provided  by  an 
absorber  covered  arm  from  the  top  or  bottom  of  the  UGAV.  A  view  of  the  400  foot  range 
at  Newport  showing  the  UGAV  and  the  ground  array  is  shown  in  Fig.  12.  It  is  seen  that 
the  range  is  relatively  free  of  scattering,  so  that  the  MIMO  performance  should  be  similar 
to  free  space  limitations. 

The  switch  at  the  ground  array  cycles  through  the  8  elements  in  a  20  millisecond  period. 
After  each  of  the  ground  array  cycles,  the  switch  at  the  UGAV  connects  to  a  new  antenna 
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((c))  Ground  Array 


Figure  11:  Element  Ground  Array 


Figure  12:  UCAV  and  Ground  Array  on  the  400  ft.  Newport  Range. 

port  of  the  16  UGAV  antenna  ports.  As  this  switching  is  completed  the  azimuth  turntable 
moves  to  a  new  position  for  a  given  pattern  cut.  At  each  position  of  the  cut,  the  complex 
transmission  coefficient  from  each  antenna  element  on  the  ground  array  to  each  element  port 
on  the  UGAV  is  recorded.  The  frequency  and  pitch  and  yaw  positions  can  be  modihed  for 
new  cuts. 

Challenges 

Fig.  10(a)  and  10(b)  illustrate  the  unique  challenges  of  airborne  MIMO.  Each  antenna  el¬ 
ement  has  only  limited  visibility  for  communication,  depending  on  the  orientations  of  the 
transceiver  pair.  This  is  further  exacerbated  by  the  lack  of  channel  state  information  at  the 
transmitter  -  RPAs  are  not  only  of  high  mobility  but  they  may  maneuver  in-flight  which 
makes  it  infeasible  to  constantly  feedback  channel  state  to  the  transmitter. 

Nevertheless,  the  channel  matrices  indicate  that  there  exists  signihcant  theoretical  through- 
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put  improvement  through  MIMO  communications  if  designed  properly.  In  the  following,  we 
describe  a  simple  variable  rate  MIMO  scheme  to  address  the  above  mentioned  challenges. 
We  demonstrate  through  numerical  simulation  that  the  proposed  scheme  is  capable  of  over¬ 
coming  the  above  challenges  and  realizing  signihcant  performance  gain. 

Measurement  Data 

We  obtain  measurement  data  for  the  following  two  conhgurations. 

LRPRM  This  conhguration  consists  of  setting  the  80  inch  wide  ground  array  at  2  feet 
above  the  ground  level.  The  polarization  of  the  antenna  elements  is  alternated  be¬ 
tween  vertical  and  horizontal  polarization.  The  UCAV  has  the  same  conhguration  as 
mentioned  in  the  previous  section. 

LRPRV  Same  as  LRPRM  except  all  the  antenna  elements  transmit  with  vertical  polariza¬ 
tion.  This  conhguration  consists  of  setting  the  80  inch  wide  ground  array  at  2  feet  above 
the  ground  level,  but  all  the  antenna  elements  transmit  with  vertical  polarization. 

4.2.2  Variable  Rate  MIMO 

While  the  channel  matrix  measurements  exhibit  potential  for  signihcant  throughput  im¬ 
provement  via  MIMO  communications,  signihcant  challenges  exist  that  need  to  be  overcome 
in  order  to  realize  this  potential.  In  particular,  the  lack  of  scattering  makes  the  channel  sus¬ 
ceptible  to  ill-conditioning  when  antenna  elements  may  be  out  of  sight  to  the  communicating 
party.  In  the  absence  of  complete  CSI  at  the  transmitter,  there  is  a  need  to  ensure  that  any 
independent  data  stream  is  not  stuck  with  an  antenna  element  that  is  out-of-sight.  This 
makes  D-BLAST  a  very  natural  candidate  for  such  applications.  D-BLAST  is  an  outage 
optimal  transceiver  architecture  for  MIMO  communication  system  [2,5].  The  architecture 
essentially  involves  independent  coding  and  modulation  of  N  data  streams  and  rotating  each 
stream  through  all  the  transmit  antennas.  Here,  N  is  the  number  of  transmit  antennas.  The 
rotation  is  performed  in  a  way  that  makes  sure  that  each  stream  experiences  the  channel 
from  all  the  transmit  antennas  and  all  levels  of  interference  i.e.,  no  interference  to  experi¬ 
encing  interference  from  all  of  the  other  A^  —  1  streams  at  the  receiver  as  each  stream  rotates 
in  the  space-time  domain. 

The  verbatim  application  of  D-BLAST,  however,  is  still  inadequate  when  channel  ma¬ 
trices  are  highly  ill-conditioned.  The  existence  of  antennas  that  experience  channel  outage 
due  to  blockage  will  drag  down  the  overall  performance.  From  a  theoretical  viewpoint,  the 
independent  data  streams  accommodated  in  a  MIMO  channel  should  be  no  larger  than  the 
effective  rank  of  the  matrix.  The  challenge  is,  in  the  absence  of  complete  CSI  at  the  trans¬ 
mitter,  how  to  utilize  the  advantage  of  D-BLAST  while  being  able  to  handle  the  potential 
rank  dehciency  of  channel  matrices. 

We  propose  a  simple  scheme  of  variable  rate  MIMO  -  the  variable  rate  is  not  only  achieved 
though  the  traditional  means  of  controlling  the  coding  rate  as  well  as  modulation  order,  but 
the  number  of  independent  data  streams  is  also  adapted  through  per  antenna  spreading. 
The  modihed  D-BLAST  architecture  is  illustrated  in  Fig.  5.  Without  loss  of  generality,  we 
assume  the  MIMO  system  has  N  transmit  and  N  receive  antennas.  At  the  transmitter  side. 
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spreading  is  done  on  a  per  antenna  basis  with  a  length  M  Walsh  code  for  each  antenna 
where  M  <  N.  As  such,  a  group  of  N/M  antennas  will  share  the  same  Walsh  code  and  the 
effective  number  of  data  streams  is  reduced  to  N/M.  At  the  receiver  side,  de-spreading  is 
embedded  in  the  D-BLAST  receiver  to  facilitate  successive  interference  cancellation. 

As  an  illustration,  consider  an  8  x  8  MIMO  system.  Let  ao,ai,...,a7  denote  the  8 
transmit  antennas.  Let  each  symbol  from  any  4  antennas,  say  00,01,02,03  be  spread  by 
[1  1]^.  Similarly,  let  each  symbols  from  the  remaining  4  antennas,  04,05,05,07  be  spread 
using  [—1  1]^.  Notice  that  [1  1]^  and  [—1  1]^  are  orthogonal  to  each  other.  At  the  receiver, 
de-spreading  effectively  reduces  the  number  of  mutually  interfering  streams  from  8  to  4. 

The  illustration  in  the  previous  paragraph  implicitly  assumes  that  the  transmitter  and 
the  receiver  have  agreed  on  some  spreading  factor,  M,  and  the  corresponding  Walsh  codes. 
As  the  transmitter  does  not  possess  any  CSI,  our  assumption  holds  true  only  if  there  exists  a 
very  low  rate  feedback  link  from  the  receiver  to  the  transmitter.  The  receiver,  which  knows 
the  CSI,  can  feedback  the  spreading  factor  to  the  transmitter,  which  then  can  adapt  the 
transmission  strategy  to  the  channel  conditions.  For  example,  in  the  8x8  case,  one  needs 
a  total  of  2  bits  for  the  spreading  factor  of  1  (no  spreading),  2,  4,  and  8  where  the  last  one 
corresponds  to  essentially  a  beamforming  strategy. 

In  the  next  section  we  simulate  the  variable  rate  MIMO  scheme  over  measured  channels 
to  examine  the  performance  of  the  proposed  variable  rate  MIMO  D-BLAST  for  low  rank 
channels. 


((b))  Receiver 

Figure  13:  Variable  rate  MIMO  transceiver  block  diagram 


4.3  Results  and  Discussion 

We  assess  the  performance  of  the  variable  rate  D-BLAST  scheme  by  comparing  it  with  the 
original  D-BLAST  for  the  following  two  channel  measurement  conhgurations. 
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•  All  Vertical  Polarization  (VP). 

•  Mixed  Polarization  (MP). 

As  mentioned  earlier,  we  are  motivated  to  use  VP  and  MP  to  see  the  effect  of  polarization 
mixing  on  MIMO  capacity.  The  channels  obtained  using  VP  tend  be  rank  dehcient.  On  the 
other  hand,  channels  obtained  using  MP  tend  to  have  a  higher  rank  than  that  of  VP.  This  is 
because  the  two  polarizations  effectively  serve  to  create  two  orthogonal  channels.  This  leads 
to  more  well-conditioned  MP  channel  matrices  than  the  corresponding  VP  ones. 

For  each  of  the  above  two  conhgurations  we  have  obtained  about  13,000  channel  measure¬ 
ments,  each  of  which  corresponds  to  a  unique  combination  of  azimuth  and  elevation.  During 
simulation,  we  select  10%  of  the  channels  at  random,  and  measure  the  average  outage  over 
these  channels. 


Figure  14:  Outage  rate  for  all  vertical  polarization  at  rate  4  bits/s/Hz. 


Figure  15:  Outage  rate  for  mixed  polarization  at  rate  4  bits/s/Hz. 

Fig.  14  compares  the  outage  behavior  of  the  variable  rate  D-BLAST  scheme  with  the 
original  D-BLAST  at  a  hxed  rate  of  4  bits/s/Hz.  The  variable  rate  D-BLAST  scheme  uses 
a  spreading  factor  of  2,  QPSK  modulation  and  1/2  rate  LDPC  code.  While  the  D-BLAST 
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scheme  uses  BPSK  modulation  and  1/2  rate  LDPC  code.  Clearly,  at  low  outage  probabilities, 
the  variable  rate  D-BLAST  scheme  provides  an  SNR  gain  of  about  4  dB  over  the  original 
D-BLAST. 

Fig.  15,  on  the  other  hand,  shows  that  the  variable  rate  D-BLAST  scheme  is  inefficient  for 
mixed  polarization  channels.  The  original  D-BLAST  outperforms  variable  rate  D-BLAST 
by  about  6  dB  SNR  at  low  outage  probabilities. 

The  above  simulation  results  show  that  for  the  low  rank  channel  matrices  obtained  by 
only  vertical  polarization,  the  variable  rate  D-BLAST  scheme  outperforms  D-BLAST,  while 
for  well  conditioned  channel  matrices,  the  original  D-BLAST  architecture  suffices. 

4.4  Conclusion 

In  this  chapter  we  have  introduced  a  new  scheme  to  transmit  data  over  MIMO  channels 
which  may  be  ill-conditioned.  In  the  absence  of  channel  state  information  at  the  transmitter 
the  scheme  overcomes  this  challenge  by  per  antenna  spreading  in  the  D-BLAST  architecture, 
thereby  effectively  reducing  the  number  of  independent  data  streams.  Using  real  measure¬ 
ment  channels  we  demonstrated  through  simulations  that  one  can  achieve  signihcant  SNR 
gain  over  the  original  D-BLAST  for  rank  dehcient  channel  matrices. 
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GnuRadio  MIMO  Implementation 
Using  USRP  N210 


5.1  Introduction 

To  facilitate  experimental  validation  of  variable  rate  MIMO  we  have  developed  a  software 
defined  2x2  MIMO  system  using  GNU  Radio/USRP  platform  [1,13]-  In  this  setup  the  GNU 
Radio  based  software  performs  all  the  baseband  signal  processing,  whereas  USRPs  act  as 
the  wireless  interfaces.  To  form  a  2  antenna  terminal,  two  USRP  N210  are  chained  together 
using  a  vendor  supplied  cable.  Thereby,  a  total  of  4  USRPs  combined  with  the  GNU  Radio 
based  baseband  signal  processing  constitute  our  2x2  MIMO  implementation.  Nevertheless, 
the  software  implementation  is  scalable  to  higher  antenna  dimensions  under  the  constraint 
that  the  number  of  transmit  antennas  is  equal  to  that  of  receive  antennas  and  some  changes 
to  the  symbol  timing  estimator  in  the  receive  chain. 

An  important  feature  of  this  implementation  is  the  decoupling  of  GNU  Radio  with  the 
DBLAST  encoding/decoding  functionality:  The  DBLAST  decoder/encoder  C++  library  is 
developed  separately  which  is  then  used  as  a  shared  library  by  the  GNU  Radio  modules. 
Such  separation  increases  the  potential  of  code  re-usability.  We  will  revisit  this  aspect  of  our 
implementation  in  Section  5.2.2.  In  the  next  section  we  provide  high  level  block  diagrams 
that  represent  the  implementation  approach. 

5.2  Methods,  Assumptions  and  Procedures 

5.2.1  System  Overview 

Fig.  16  shows  the  transmitter  block  diagram.  As  shown,  an  input  byte  stream  is  first  de¬ 
multiplexed  and  then  given  to  the  DBLAST  Encoder.  The  DBLAST  encoder  performs 
channel  encoding,  digital  constellation  mapping  and  the  diagonal  layering  of  the  two  incom¬ 
ing  streams.  The  two  output  streams  of  the  encoder  are  then  fed  to  two  upsampling  RRC 
Filters  (root  raised  cosine).  The  outputs  of  these  Liters  are  then  fed  to  two  separate  USRP 
N210s,  which  share  a  common  10  MHz  reference  clock. 

Fig.  17  shows  the  structure  of  the  receiver.  The  outputs  from  the  two  synchronized 
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Figure  16:  Transmitter  Block  Diagram 


USRP  N210s  are  RRC  filtered  and  processed  to  obtain  the  symbol  timing  and  then  sampled 
accordingly.  The  sampled  output  is  then  compensated  for  the  carrier  frequency  offset  (CFO) 
between  the  transmit  and  receive  USRPs.  The  CFO  corrected  output  is  then  given  to  the 
DBLAST  Decoder  which  performs  MMSE  SIC  detection  on  the  received  signal. 

5.2.2  CH — h  MIMO  Library 

As  mentioned  previously,  we  have  implemented  DBLAST  encoding  /  decoding  as  a  separate 
C++  library  which  is  then  used  by  our  GNU  Radio  implementation.  The  library  can  be  used 
independently  of  GNU  Radio  and  makes  use  of  Eigen  library  [14]  for  matrix  manipulations. 
The  DBLAST  decoder  therein  implements  MMSE  SIG  algorithm.  The  library  along  with 
its  documentation  has  been  delivered  to  AFRL. 


5.3  Results  and  Discussion 

Using  the  GNU  Radio/USRP  software  radio  platform  we  successfuly  implemented  a  DBLAST 
transciever.  We  also  performed  over  the  air  tests  with  success  and  thereby  validated  the  im¬ 
plementation.  An  independent  C++  library  that  provides  DBLAST  encoder/decoder  func¬ 
tionality  was  developed  and  delivered  as  well. 
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Conclusion 


In  this  report  we  studied  problems  relating  to  airborne  MIMO  communication  and  delivered 
relevant  MATLAB,  GNU  Radio  and  C++  based  simulation  and  experimentation  software. 

Firstly,  we  studied  the  problem  of  channel  tracking  in  dynamic  conditions  entailed  in 
airborne  MIMO  communication.  In  doing  so,  we  developed  a  scheme  to  estimate  and  keep 
the  channel  state  information  current  using  payload  data.  We  analyzed  the  scheme  and 
provided  a  closed  form  expression  to  determine  optimal  symbol  length  needed  to  estimate 
the  channel  in  a  continuously  fading  Rayleigh  environment.  In  addition,  we  have  provided 
a  channel  tracking  simulator,  developed  in  MATLAB,  to  quickly  test  different  scenarios. 

Secondly,  by  studying  the  real  measurement  data,  we  found  out  that  the  wireless  channel 
matrices  offered  therein  are  severely  rank  dehcient.  Traditional  MIMO  architectures,  which 
are  developed  to  take  advantage  of  well  conditioned  MIMO  channel  matrices,  fail  to  take 
advantage  of  such  rank  dehcient  channels.  We  developed  a  Variable  Rate  MIMO  scheme  to 
overcome  the  problem  of  communication  over  rank  dehcient  channels.  Through  simulations 
we  showed  that  the  developed  scheme  signihcantly  outperforms  D-BLAST  for  rank  dehcient 
channel  matrices. 

Finally,  we  provided  a  GNU  Radio  based  2x2  D-BLAST  implementation,  to  facilitate 
experimentation.  Further,  we  developed  a  G++  library  that  implements  D-BLAST  encoder 
and  decoder. 
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Symbols,  Abbreviations  and 
Acronyms 


0  zero  vector  with  appropriate  dimension  depending  on  the  context 
Ifc  identity  matrix  of  dimension  with  k  columns  and  k  rows. 

X  matrix  denoted  by  uppercase  bold  letter 
X  column  vector  denoted  by  lowercase  bold  letter 
*  complex  conjugate  operator 
+  Hermitian  matrix  operator 
^  transpose  operator 

r  X  t  wireless  communication  system  with  r  receive  antennas  and  t  transmit  antennas 
X  scalar  denoted  by  lower  case  letter 

BPSK  binary  phase  shift  keying 

CSI  channel  state  information 

D-BLAST  diagonal  Bell  laboratories  layered  space-time 
GUI  graphical  user  interface 
LDPC  low  density  parity  check 

MIMO  multiple  input  multiple  output 
MSE  mean  square  error 

QAM  quadrature  amplitude  modulation 
QPSK  quadrature  phase  shift  keying 

RF  radio  frequency 
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RPA  remotely  piloted  aircraft 


SNR  signal  to  noise  power  ratio 
UCAV  unmanned  combat  air  vehicle 
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Appendix  A 

Proof  of  Theorem  1 


Proof.  Using  the  assumptions  mentioned  in  the  text,  we  first  evaluate  ^(ApAU. 


^(ApA! 


=  ^ 


j2'Kfk{N-p)Ts  _  j2iTfkNTs 


.k=l 


-j2'nfi{N-q)Ts  _  -j2nfiNTs 


1=1 


= 


E  lAP  (< 


J2wfk{N-p)Ts  _  j2TTff,NTs 


.k=l 


+<P 


^-3‘^'^fk{N-q)Ts  _  -j2-KfkNTs 


k^l 


(^Q-p'^fl{N-q)Ts  _  ^-j2-KfiNTs 


s[(l  -  (l-e 


,j2-nfuqTs 


E 

.k=l 


=  s 


(l  - 


oP'^fkqTa 


=  1  -  smc{2TT fopTs)  -  smc{2TTfDqTs) 

+smc{27rf d{p  -  g)Tj,  (17) 

where  sinc(a;)  =  Third  equality  is  because  the  second  expectation  in  the  second 

step  is  zero. 

From  the  above  expression  we  see  that, 

^(ApAp  =  ^(A,A;).  (18) 

Also,  putting  p  =  q  in  expression  (17),  we  get, 

^{ApA*)  =  2(1  -  sinc(27r/DpTj).  (19) 
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We  now  simplify  the  expectation  term  in  equation  (6).  Using  brute  force  we  can  show 
that, 

‘•[(S*-)  (£«)] 

■  2  2 

=  2  smc{27r fokTs)  ,  (20) 

.k=l  k=l 

bl:,*-)  (5«)1 

■  3  3 

=  2  '^k  -'^k  smc{2n fokTs)  .  (21) 

.fc=i  k=i 

We  now  simplify  the  expectation  term  in  equation  (6).  We  use  proof  by  induction. 
Suppose  the  following  is  true, 

'•[(£*■)  (sbl 

r  N  N  1 


Now, 


=  2  smc{27r fokTs) 

.k=l  k=l 

niseis 

=  ^  f X!  + A) 


Let,  (Ef=iAfc)  =D. 


- /N+l  \  /W+1 
.\fc=i  /  \fc=i 


—  S'  [D  +  Ajv+i)  (^D*  + 

=  S{DD*)  +  S{DA%^,)  +  S{A^+,D*) 

+S{A]\[^iA*]^_^_i). 

We  know  S{DD*)  from  equation  (22).  We  need  to  hud  the  last  three  terms. 
S{DA*^^,) 

=  S  (Ai  +  A2  +  . . .  +  A7v)A^_,_;^) 

=  ^(AlA^^;^)  +  ^(A2A^_,_;^)  +  .  .  .  +  ^(AlA^_,_;^) 

N  N 

=  A  -  y])  sinc(27r/D/cTj  -  J])  sinc(27r/D(A  +  1)T^) 

N 

+  y])  sinc(27r/D(A  +  1  -  k)Ts) 


=  N-N  smc{27rf d{N  +  1)T,). 
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From  equation  (19), 


=  2(1  -  sinc(27r/D(A  +  1)TJ). 

From  equations  (18)  and  (27), 

^(DA^+J  =  ^(A^+iD*). 

From  equations  (22),  (25),  (29),  (30)  and  (31),  we  get 


(30) 


(31) 


/N+l  \  /N+1  \ 

i:  aO  x  Ai 

\A:=1  /  \A:=1  V 


=  2 


■  N  N 

E^-E  k  smc{27r fnkTs 

.k=l  k=l 

+2{N  -  Asinc(27r/D(iV  +  1)T,)) 

+2(1  -  sinc(27r/D(iV  +  l)Ts)) 
W+l  Af+l 

^  smc{27r fokTs) 


=  2 


k=l  k=l 


(32) 

(33) 


The  last  term  in  the  above  expression  is  of  the  form  J2k=i  ksmc{ck),  where  c  is  a  constant. 


Therefore, 


ck 


N 

k  sinc(ck) 

k=l 

^  sinc/c 

=  2^k 

k=l 
1  ^ 

=  -  sin  ck 

^  k=i 

1  .  fNc\  . 
=  -  sm  —  sm 

c  \  2  J 


(A  +  l)c\  /c 

2  J  V2 


(34) 


N  \  /  ^  '' 

k=l  /  \k=l  / 

■  N  N 

k  —  k  sinc{2'K fokTg 

—  1  k^—  1 

sm{7r foNTs)  sm{7ifD{N  +  1)TJ 


=  2 

=  N{N  +  1)  - 


tt/dTs  sm{Tr foTs) 


(35) 

(36) 
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Chapter  1 


Class  Index 

1.1  Class  Hierarchy 

This  inheritance  list  is  sorted  roughly,  but  not  completely,  alphabetically: 

mimo::dblast::base  .  5 

mimo::dblast::decoder .  6 

mimo::dblast::encoder .  12 

mimo::enc_dec  .  11 

mimo::decoder .  8 

mimo::encoder .  14 

mimo::mod_demod .  14 

mimo::demodulator .  9 

mimo::demodulator_bpsk .  10 

mimo::demodulator_qpsk .  11 

mimo::modulator .  15 

mimo::modulator_bpsk  .  16 

mimo::modulator_qpsk  .  17 

mimo::repacker .  18 
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Class  Index 
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Chapter  2 


Class  Index 

2.1  Class  List 

Here  are  the  classes,  streets,  unions  and  interfaces  with  brief  descriptions: 

mimo::dblast::base  .  5 

mimo::dblast::decoder  .  6 

mimo::decoder  .  8 

mimo::demodulator .  9 

mimo::demodulator_bpsk .  10 

mimo:;demodulator_qpsk .  11 

mimo::enc_dec  .  11 

mimo::dblast::encoder  .  12 

mimo::encoder  .  14 

mimo::mod_demod .  14 

mimo::modulator  .  15 

mimo::modulator_bpsk  .  16 

mimo::modulator_qpsk  .  17 

mimo::repacker .  18 
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Chapter  3 


Class  Documentation 


3.1  mimo;;dblast::base  Class  Reference 

Inheritance  diagram  for  mimo::dblast::base: 


Public  Member  Functions 

•  unsigned  ntx  ()  const 

•  unsigned  nrx  ()  const 

•  unsigned  message_length_bits  () 

•  unsigned  block_length_bits  () 

•  unsigned  block_length_symbs  () 

•  unsigned  bits_per_symbol  () 

•  unsigned  sub_block_length_symb  ()  const 

•  unsigned  num_sub_blocks  (const  unsigned  news) 


base  -  Base  Class  for  DBLAST  encoder  and  decoder 
Parameters 


ntx 

-  Number  of  transmit  antennas 

nrx 

-  Number  of  Receive  antennas 

mod 

-  Moduiator  object 

enc 

-  Encoder  object 
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Class  Documentation 


Returns 

base  object 

•  base  (const  unsigned  ntx,  const  unsigned  nrx,  mimo::modulator  &mod,  mimo- 
::encoder  &enc) 

Protected  Member  Functions 

•  void  set_num_sub_block  () 

Protected  Attributes 

•  const  unsigned  d_ntx 

•  const  unsigned  d_nrx 

•  unsigned  d_sbl 

•  unsigned  d_bls 

•  mimo::repacker  *  d_packer 

•  mimo::modulator  &  d_mod 

•  mimo::encoder  &  d_enc 

^base  -  defauit  destructor 

Returns 

void 


•  ~base  {) 

•  void  check_args  () 

•  void  set_sbl  () 

The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/dblast.h 

•  /home/kapil/libmimo/lib/dblast.cc 


3.2  mimo;;dblast;;decoder  Class  Reference 

Inheritance  diagram  for  mimo::dblast::decoder: 
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3.2  mimo::dblast::decoder  Class  Reference 
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Public  Member  Functions 

decoder  -  DBLAST  Decoder 


Parameters 


ntx 

-  Number  of  transmit  antennas 

nrx 

-  Number  of  receive  antennas 

mod 

-  Modutator  object 

demod 

-  Demoduiator  object 

enc 

-  Encoder  Object 

dene 

-  Decoder  Object 

ninput_cws 

-  Number  of  input  codewords  in  a  block 

Returns 

decoder  object 

•  decoder  (const  unsigned  ntx,  const  unsigned  nrx,  mimo::modulator  &mod, 
mimo::demodulator  Sdemod,  mimo::encoder  &enc,  mimo::decoder  &denc, 
const  unsigned  ninput_cws) 


set_ninput_cws  -  Allocates  memory  for  use  during  encoding. 


Parameters 


ninput_cws 


-  Number  of  input  codewords 


Returns 

void 

•  void  set_nlnput_cws  (const  unsigned  news) 

decode  -  Decodes  received  DBLAST  signal 


Parameters 


inx 

-  Matrix  of  dimenstion  nsymb  x  ntx  where  nsymb  is  the  number  of  sym¬ 
bols  time  slots  and  ntx  is  the  number  of  transmit  antennas 

Hhat 

-  ntx  nrx  channel  matrix 

snr 

-  Signai  to  noise  ratio. 

Returns 

void 

•  void  decode  (const  Eigen::MatrixXcf  &inx,  const  Eigen:;MatrixXcf  &Hhat, 
const  float  &snr) 

output_mat  -  Get  decoded  bits 
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Returns 

Returns  constant  reference  to  Matrix  of  type  mimocMatrixXuoharb.  Each  col¬ 
umn  corresponds  to  a  deooded  codeword. 

•  const  MatrixXuchar  &  output_mat  {)  const 


output  -  Returns  the  output  generated  by  the  decode  function 


Parameters 


out 


Output  veotor  to  filled  with  decoded  bits 


Returns 

void 

•  void  output  {vector<  uchar  >  &out)  const 

~decoder  -  default  destructor 

Returns 

void 

•  ^decoder  () 

The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/dblast.h 

•  /home/kapil/libmimo/lib/dblast.cc 

3.3  mimo::decoder  Class  Reference 

Inheritance  diagram  for  mimo::decoder: 


Public  Member  Functions 

encoder  -  Generic  decoder  class 
Parameters 


n 

-  Codeword  block  length  in  bits 

k 

-  Codeword  message  length  in  bits 
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3.4  mimo::demodulator  Class  Reference 
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Returns 

decoder  object 

•  decoder  (unsigned  n,  unsigned  k) 

•  ^decoder  () 

•  virtual  void  decode  (const  vector<  uchar  >  &x,  vector<  uchar  >  &y) 

•  virtual  void  decode  (const  uchar  *inp,  const  unsigned  inlen,  uchar  *outp,  const 
unsigned  outlen) 


The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/enc_dec.h 

•  /home/kapil/libmimo/lib/enc_dec.cc 


3.4  niimo::demodulator  Class  Reference 

Inheritance  diagram  for  mimo::demodulator: 


Public  Member  Functions 

~modulator_qpsk  -  default  destructor 

Returns 

void 

•  demodulator  () 

•  demodulator  (const  vector<  complex<  float  >  >  &constellation) 

•  ^demodulator  () 

•  virtual  unsigned  decode  (const  complex<  float  >  Ssample) 

•  virtual  void  decode  (const  complex<  float  >  ^sample,  uchar  *out,  const  un¬ 
signed  vien) 

•  virtual  void  decode  (const  complex<  float  >  ^sample,  unsigned  *out,  const 
unsigned  vIen) 

•  virtual  void  decode  (const  vector<  complex<  float  >  >  Ssamples,  vector< 
unsigned  >  &out) 
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Protected  Attributes 

•  vector<  float  >  d_distances 

•  unsigned(demodulator::>K  decode_sample_ptr  )(const  complex<  float  >  &) 
The  dooumentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/inelude/mod_demod.h 

•  /home/kapil/libmimo/lib/mod_demod.cc 

3.5  mimo::demodulator_bpsk  Class  Reference 

Inheritance  diagram  for  mimo::demodulator_bpsk: 


Public  Member  Functions 

~modulator_qpsk  -  default  destructor 

Returns 

void 

•  demodulator_bpsk  () 

•  ~demodulator_bpsk  () 

•  unsigned  decode  (const  complex<  float  >  Ssample) 

•  void  decode  (const  complex<  float  >  ^sample,  uchar  *out,  const  unsigned 
vien) 

•  void  decode  (const  complex<  float  >  ^sample,  unsigned  *out,  const  unsigned 
vIen) 

•  void  decode  (const  vector<  complex<  float  >  >  Ssamples,  vector<  unsigned 
>  &out) 


The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/mod_demod.h 

•  /home/kapil/libmimo/lib/mod_demod.cc 
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3.6  mimo::demodulator_qpsk  Class  Reference 
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3.6  niimo::demodulator  qpsk  Class  Reference 

Inheritance  diagram  for  mimo::demodulator_qpsk: 


Public  Member  Functions 

~modulator_qpsk  -  default  destructor 

Returns 

void 

•  demodulator_qpsk  () 

•  ~demodulator_qpsk  () 

•  unsigned  decode  (const  complex<  float  >  Ssample) 

•  void  decode  (const  complex<  float  >  ^sample,  uchar  *out,  const  unsigned 
vien) 

•  void  decode  (const  complex<  float  >  ^sample,  unsigned  *out,  const  unsigned 
vien) 

•  void  decode  (const  vector<  complex<  float  >  >  &samples,  vector<  unsigned 
>  &out) 


The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/mod_demod.h 

•  /home/kapil/libmimo/lib/mod_demod.cc 


3.7  mimo::enc_dec  Class  Reference 


Inheritance  diagram  for  mimo::enc_dec: 
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Public  Member  Functions 

•  unsigned  block_length_bits  () 

•  unsigned  message_length_bits  () 


enc_dec  -  Base  class  for  encoder  and  decoder  class 
Parameters 


n 

-  Codeword  block  length  in  bits 

k 

-  Codeword  message  length  in  bits 

Returns 

enc_dec  object 

•  enc_dec  (unsigned  n,  unsigned  k) 

•  ~enc_dec  () 


Protected  Member  Functions 

•  string  int_to_string  (int  x) 


Protected  Attributes 

•  unsigned  d_n 

•  unsigned  d_k 


The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/enc_dec.h 

•  /home/kapil/libmimo/lib/enc_dec.cc 


3.8  mimo::dblast::encoder  Class  Reference 


Inheritance  diagram  for  mimo::dblast::encoder: 
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Public  Member  Functions 


encoder  -  DBLAST  Encoder  Constructor 
Parameters 


ntx 

-  Number  of  transmit  antennas 

nrx 

-  Number  of  receive  antennas 

mod 

-  Reference  to  modulator  object 

enc 

-  Reference  to  encoder  object 

ninput_cws 

-  Number  of  Input  codewords  In  a  block 

Returns 

DBLAST  Encoder  Object 

•  encoder  (const  unsigned  ntx,  const  unsigned  nrx,  mimo::modulator  &mod, 
mimo::encoder  &enc,  const  unsigned  ninput_cws) 

^encoder  -  default  destructor 
Returns 

void 

•  ^encoder  () 


set_nlnput_cws  -  Allocates  memory  for  use  during  encoding. 


Parameters 


ninput_cws 


-  Number  of  input  codewords 


Returns 

void 

•  void  set_nlnput_cws  (const  unsigned  ninput_cws) 


encode  -  Encodes  the  input  stream  of  bits  for  DBLAST  transmission 


Parameters 


inx 


-  Input  vector  of  type  mimo::uchar.  Veotor  length  must  be  multiple  of  the 
codeword  block  length  In  bits. 


Returns 

void  (To  aocess  the  result,  call  output  function.) 
•  void  encode  (const  vector<  uchar  >  &inx) 


output 
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Returns 

returns  reference  to  the  output  matrix  of  type  Eigen::MatrixXcf 

•  const  Eigen  ::MatrixXcf  &  output  () 

The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/dblast.h 

•  /home/kapil/libmimo/lib/dblast.cc 


3.9  mimo::encoder  Class  Reference 

Inheritance  diagram  for  mimo::encoder: 


Public  Member  Functions 


encoder  -  Generic  encoder  class 
Parameters 


n 

-  Codeword  block  length  In  bits 

k 

-  Codeword  message  length  In  bits 

Returns 

encoder  object 

•  encoder  (unsigned  n,  unsigned  k) 

•  ^encoder  () 

•  virtual  void  encode  (const  vector<  uchar  >  &x,  vector<  uchar  >  &y) 

•  virtual  void  encode  (const  uchar  *inp,  const  unsigned  inlen,  uchar  *outp,  const 
unsigned  outlen) 

The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/enc_dec.h 

•  /home/kapil/libmimo/lib/enc_dec.cc 


3.10  niimo::mod  demod  Class  Reference 

Inheritance  diagram  for  mimo::mod_demod: 
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Public  Member  Functions 

•  mod_demod  (const  vector<  complex<  float  >  >  Sconstellation) 

•  int  mod_order  () 

•  unsigned  bits_per_symbol  {) 

•  virtual  vector<  complex<  float  >  >  points  {) 

Protected  Member  Functions 

•  void  set_bps  () 

Protected  Attributes 

•  int  d_mod_order 

•  unsigned  d_bps 

•  vector<  complex<  float  >  >  d_points 

The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/mod_demod.h 

•  /home/kapil/libmimo/lib/mod_demod.cc 


3.11  mimo::modulator  Class  Reference 

Inheritance  diagram  for  mimo::modulator: 
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Public  Member  Functions 

•  modulator  (const  vector<  complex<  float  >  >  Sconstellation) 

•  virtual  complex<  float  >  encode  (const  unsigned  idx) 

•  virtual  void  encode  (const  unsigned  idx,  complex<  float  >  *out) 

•  virtual  void  encode  (const  unsigned  *idx,  complex<  float  >  *out,  const  unsigned 
vien) 

•  virtual  void  encode  (const  mimo::uchar  *idx,  complex<  float  >  *out,  const  un¬ 
signed  vIen) 

•  virtual  void  encode  (const  vector<  unsigned  >  &idx,  vector<  complex<  float  > 
>  &out) 

The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/mod_demod.h 

•  /home/kapil/libmimo/lib/mod_demod.cc 


3.12  mimo;;modulator_bpsk  Class  Reference 

Inheritance  diagram  for  mimo::modulator_bpsk: 


Public  Member  Functions 

modulator_bpsk  -  BPSK  modulator  class 
Returns 

modulator_bpsk  object 

•  modulator_bpsk  () 

~modulator_bpsk  -  default  destructor 
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Returns 

void 

•  ~modulator_bpsk  {) 

The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/mod_demod.h 

•  /home/kapil/libmimo/lib/mod_demod.cc 


3.13  mimoiimodulator  qpsk  Class  Reference 

Inheritance  diagram  for  mimo::modulator_qpsk: 


Public  Member  Functions 

modulator_qpsk  -  Grey  coded  QPSK  modulator  class 
Returns 

modulator_qpsk  object 

•  modulator_qpsk  () 

~modulator_qpsk  -  default  destructor 

Returns 

void 

•  ~modulator_qpsk  {) 

The  documentation  for  this  class  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/mod_demod.h 

•  /home/kapil/libmimo/lib/mod_demod.cc 
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3.14  mimo::repacker  Class  Reference 

Public  Member  Functions 

•  repacker  (const  unsigned  k,  const  unsigned  i) 

•  void  operatorQ  (const  uchar  *inp,  const  unsigned  inlen,  uchar  *outp,  const  un¬ 
signed  outlen) 

The  documentation  for  this  ciass  was  generated  from  the  following  files: 

•  /home/kapil/libmimo/include/repacker.h 

•  /home/kapil/libmimo/lib/repacker.cc 
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